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We study the effect of external driving on the two level systems (TLSs) assumed to be a major 
obstacle in increasing the coherence time of solid state Josephson-junction qubits. We find, by use 
of a Bloch-Redfield approach, that external driving has two major effects on the TLS. The first is 
increased fluctuations between the two states of the TLS, the significance of this effect compared 
to thermal fluctuations depend on the energy splitting of the TLS compared to temperature. The 
second effect is a reduction in the intensity of the noise spectrum at low frequencies, and at the 
same time an increase in intensity around the renormalized Rabi frequency of the TLS, the driving 
frequency and at beatings between these two frequencies. Finally we study the ensemble averaged 
noise spectrum for a typical distribution of TLSs known to give origin to oc 1// noise. We find that 
strong driving leads to reduced noise at low frequencies, and therefore to an increased dephasing 
time of the qubit. However this effect is exponentially suppressed when the driving frequency is 
large compared to temperature, as we typically find for Josephson qubits. We suggest that external 
driving at frequencies much lower than the qubit frequency might be used in order to enhance the 
the qubit coherence time. 

PACS numbers: 42.50. Lc, 03.67.Lx, 03.65.Yz, 74.78.-w, 85.25.Cp 



I. INTRODUCTION 

The most fundamental problem that has to be over- 
come in order to produce a quantum computer is the iso- 
lation of its basic elements, the quantum bits (qubits), 
from its environment. Entanglement with uncontrollable 
degrees of freedom is responsible for the decay of coherent 
superpositions of qubit states. The result is irreversible 
loss of the quantum information required for operation of 
the device. Superconducting qubits based on the Joseph- 
son junction are leading candidates in the design of a 
quantum computer. They have low losses, are easily con- 
trollable by microwave pulses and can be fabricated by 
use of established integrated circuit technology. Recent 
progress in extending the decoherence time of the qubits 
has been achieved by identification of the sources of noise 
and their respective natures. This knowledge has lead to 
the development of countermeasures such as better isola- 
tion, as well as protocols to minimize the negative impact 
of the noise, see, e.g., Refs. QJ-0, or Ref. for a review. 

Bistable two level systems (TLS) existing in the tun- 
neling junction and in the amorphous substrate used to 
fabricate the qubit, are thought to be the most important 
source of decoherence in Josephson junction qubits r"~— 
see also Ref. [ID for a review. These TLSs are assumed to 
give rise to the observed oc 1// noise spectra in Joseph- 
son qubits. It is known that control and manipulation 
of Josephson qubits by use of microwave pulses unavoid- 
ably leads to driving of TLSs in the vicinity of the qubit. 
While different theories of oc 1/f noise, and their conse- 
quences have been studied in great detail, the effect of 
driving has with one exception been neglected. 

Recently, the influence of external driving on the noise 



spectra of such TLS was investigated in Ref. [TH It was 
found that the noise at low frequencies was unchanged by 
driving, while the noise at high frequencies was weakly 
reduced. In this article we calculate the noise spectra 
from a single TLS and an ensemble of driven TLSs. The 
results we obtain differ qualitatively from those obtained 
in Ref. [H|. The picture we arrive at is the following: In a 
general environment, e.g., a disordered substrate, there 
will be TLSs with a wide distribution of energy splittings 
E and relaxation rates 7. Given a driving field of fre- 
quency f2, we can divide the fluctuators into two groups. 
Those who are far from resonance and very weakly per- 
turbed by the driving field (group /), and those who are 
close to resonance with the driving field (group II). We 
find that the TLSs belonging to group 77, are strongly af- 
fected by the driving provided that the driving amplitude 
is large compared to the relaxation rate of the TLSs. The 
response to the driving can roughly be described by two 
effects. The first is saturation of the fluctuators. A two- 
level fluctuator with large energy splitting compared to 
temperature, E 3> fc^T, will in the absence of driving be 
frozen in the ground state, with a very small probability 
of switching to the exited state. By driving this fluctua- 
tor with a frequency close to resonance, the probability 
of exitation will increase and by increasing the driving 
intensity the probability for the fluctuator to be found 
in the upper state versus the lower state will eventually 
be similar, thus the TLS is saturated. A driven fluctua- 
tor will thus fluctuate (much) more rapidly between its 
upper and lower state. Therefore the noise from this fluc- 
tuator will increase. The second effect caused by driving 
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at frequencies centered around lo = and at the same 
time increased noise at higher frequencies. The driving 
results in several new peaks in the noise spectrum. Most 
pronounced are the peaks centered at the renormalizcd 
Rabi frequency, A', at the driving frequency, f2, and at 
beatings between these two frequencies. We find that the 
net effect of strong driving is a suppressed noise spectrum 
at low frequencies. 

A typical substrate used for fabrication of qubits is of- 
ten assumed to contain TLSs with a oc 1/7 distribution 
of relaxation rates, and a smooth distribution of energy 
splittings, which can be approximated as uniform.— i 12 ' 1 
This distribution is known to give rise to oc 1// noise 
at low frequencies^ For such an ensemble of fluctua- 
tors, the low frequency noise is strongly dominated by 
the fluctuators with small relaxation rates. Driving at 
high frequencies resonant with the energy splitting of the 
qubit, HQ « Eq where Eq > UbT (the energy splitting 
needs to be large compared to temperature, in order to 
avoid thermal transitions between its eigenstates) , will 
result in a significant response only from the fluctuators 
near resonance with the driving field. In the absence of 
driving, this subset of fluctuators (group II) arc frozen 
out and contribute only marginally to the ensemble aver- 
aged noise spectra, which is dominated by the fluctuators 
with small energy splittings, E < ksT. Thus suppres- 
sion of the low frequency noise from group 77 by strong 
driving only weakly influences the full ensemble-produced 
noise spectra at low frequencies, but strongly increases 
the noise at higher frequencies. 

However, our results show that while external driv- 
ing at qubit frequency for typical ensembles of fluctua- 
tors will not have significant impact on the low frequency 
noise, external low frequency driving might significantly 
suppress it . Driving at low frequencies will effect the fluc- 
tuators that contribute most strongly to the dephasing- 
producing noise felt by the qubit (i.e., those with small 
energy splitting, E). These TLSs are only weakly influ- 
enced by the saturation effect since their ratio E/ksT is 
low and correspondingly the population level in the upper 
state is already high in the absence of driving. Therefore, 
the net effect of driving on the low E fluctuators is almost 
entirely a shift in the frequency spectra from low to high 
frequencies. To us this seems like a promising method to 
reduce pure dephasing noise and thereby increase for 
the qubit. It is, however, important to note that the high 
frequency noise will be increased, specifically around the 
Rabi frequency of the driven fluctuators as well as around 
the driving frequency. One should therefore make sure 
that these frequencies lie sufficiently far from the eigen- 
frequency of the qubit in order to avoid decreasing . 
In this article we will focus on the low frequency noise, 
the noise at frequencies close to the qubit splitting need 
to be treated separately, see Refs. Il6l - fl8l . 



The rest of this article is divided in the following sec- 
tions. In Scction[TI]we will describe our model of a TLS in 
an external field and the assumptions behind it. There- 
after in Section IIIII we will derive an expression for the 
noise spectrum from a single TLS and look at different 
limiting cases. In Section IIVI we will derive an expres- 
sion for the ensemble averaged noise in the case of strong 
driving, for a particular distribution of TLS parameters 
P{E, r y) oc I/7. In Scction|V]wc will discuss the effect of 
the driven TLSs on the central qubit. Finally the results 
will be discussed in Section [VTl 



II. MODEL 

In this section we will study the dynamics of TLSs 
(fluctuators) subject to an external AC electric field, 
E ac , and a thermal environment. The nature of the two 
level systems we are interested in can, e.g., be consid- 
ered to be bistable fluctuators tunneling between dis- 
tinct charge configurations, leading to charge noise in 
the qubit. These charge fluctuators might be attributed 
to tunneling of charges between either localized impu- 
rity states, between localized impurity states and metal- 
lic electrodes, or between different charge configurations 
in a dielectric material . 10 ' 19 We model the charge config- 
urations associated with each state of a given TLS by its 
effective dipolc moment p. In order to capture the action 
of the environment (e.g., thermal phonons) responsible 
for relaxation and decoherence of the TLS, we apply the 
Bloch-Rcdfield approach^Sr— We assume that the inter- 
action between different TLSs is weak compared to the 
coupling to the thermal bath, such that eventual corre- 
lations between the TLSs are neglected. Furthermore we 
assume that the TLSs couple sufficiently weakly to the 
qubit compared to other degrees of freedom in the en- 
vironment that neglection of the qubit is justified when 
studying the dynamics of the TLS. This allows us to use 
a perturbative approach when treating the effect of the 
TLS(s) on the qubit. 



A. Hamiltonian 

Our Hamiltonian for the TLSs closely follows that of 
Ref . Il3l. A fluctuator, e.g., a particle in a double well po- 
tential with associated dipole moment p, can be modeled 
as a two level system with tunneling matrix element Ao 
and asymmetry energy A. The Hamiltonian of this TLS 
in an applied electric field, E QC is then H(t) = H +Hi (t), 
where H = \ (Ar 2 + A t x ) and H x (t) = -T z pE ac (*). 
Here t x , z are the Pauli matrices and E ac (i) = E ac cosfJi 
is an AC electric field of angular frequency f2 coupling to 
the electric dipole moment of the TLS. Furthermore, the 
TLS interacts with the qubit through Hf-q and couples 
to the environment through Hp- env . 

By diagonalization of Hq, the Hamiltonian in the en- 
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crgy cigcnbasis becomes 



H — Ho + Hi(t) + Hp-Q + Hp- env , 
H = \eo,, 



H\{t) = -r](Aa z + A <r x ) cosilt, 
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While not valid in general, this approximation is believed 
to be valid when the decoherence is isotropic^ The gen- 
eral relationship T2 < 1T\ can be derived from the master 
equation approach] 20 ' 24 Thus by making this simplifying 
assumption, asymmetry of the relaxation behavior of the 
TLS are left out. However we believe that these details 
are not of crucial importance for the results derived con- 
cerning the qubit's decoherence due to the TLSs. 



where E 
jj, z and cr : 



\J A 2 + Aq and r\ = p~E ac /E. The matrices 
are Pauli matrices acting in the eigenbasis of 
the qubit and the TLS, respectively and v is the qubit- 
fluctuator coupling parameter. The TLS-qubit coupling 
will be neglected when treating the dynamics of the TLS, 
assuming it is weak compared to other terms. But it is, 
of course, important with regards to the decoherence of 
the qubit. We note that the situation when the qubit 
and the TLS have very close splittings is an exception. 
Then the interaction is strong. See, e.g., Refs. Ilftjl8l 

Rather than specifying the explicit nature of the cou- 
pling to the environment, Hp-env, we make use of the 
Bloch-Redfield equation, where the environment enter as 
damping terms, seeking to relax the density matrix to- 
wards its thermal equilibrium value. The Bloch-Redfield 
equations for the density matrix elements of the two level 
system in the eignenbasis of Ho are^ 
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= T 2 _ can be derived from perturbation the- 
ory^ and the equilibrium density matrix, p e ^ p,(T) is in- 
troduced phcnomenologically in order to achieve relax- 
ation towards thermal equilibrium. 

By use of the rotating wave approximation, we can 
simplify the first two terms of our Hamiltonian, Eq. ([2} , 
obtaining 
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Inserted into the Bloch-Redfield equation Eq.© we find 
that the time evolution of the elements of the density ma- 
trix is governed by the following set of differential equa- 
tions: 
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III. SINGLE TLS 

In this section we will first solve the equations of mo- 
tion for a single TLS and then proceed to find its noise 
spectrum. This we will analyze later when we study the 
influence of the TLS(s) on the qubit. Using the nota- 
tions N = p++, Re / = R and Im/ = / one can cast the 
Bloch-Redfield equation, Eq. ([5]), in the form 



N = -AI - 7 (N - N cq ) , 
i = A(N - 1/2) - ~/I + zR, 
R = -zI-jR. 

Solution of Eqs. © can be written as 
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where Nq, Ro and Iq are the initial values of N, R and 
J, respectively. The solution of the homogeneous part of 
the equation is given by 
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where A' = V A 2 + z 2 is the renormalized Rabi fre- 
quency. The particular solution is given by 
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We note that by setting t — > 00 in Eq.® it is possible to 
directly read out the steady state solution. 



Here we note that p^ = pL + and p = 1 — /?++. To 

avoid the explicit time dependence we make the transfor- 
mation / = e~ int p |_ and /* = e' n *p+- • We also intro- 
duce the Rabi frequency A = r/Ao/h and the deviation 
from resonance z = E/H — Q. Furthermore we make the 
approximation for the relaxation rates 7 = 1/T\ = I/T2. 



A. Noise spectrum from a single TLS 

Given the above specified qubit-TLS coupling, Eq.©, 
the TLS is only responsible for pure dephasing of the 
qubit (T 2 processes) and cannot induce transitions be- 
tween the eigenstates of the p z operator. The dynamics 
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of the TLS leads to uncontrolled fluctuations in the en- 
ergy splitting of the qubit, leading to an uncertainty in 
its phase. Alternatively, it leads to entanglement both 
directly to the fiuctuators and indirectly to the environ- 
ment of the fiuctuators. In this article we will analyze 
the effect of the fiuctuators on the qubit through the 
two-time correlation function of the operator responsible 
for the noise in the qubit energy splitting^ We define it 
as 

G(h,t 2 ) = ([q(t 2 ) - q(t 2 )] [q(h) - q(h)}) 
= ^2{q(h) - q(t2))\ qj ( tl )-q( tl ) 

j 

x [ qj ( tl ) - q^Piqjih) - q{h)]. (10) 

Here qj(t) is a realization of a measurement at time t 
of the operator vt z , giving the variation in the qubit's 
energy splitting due to its interaction with a TLS. While 
q(t) = (q(t)) = J2j qj(t)P[<lj(t)] 15 ^e ensemble average 
of q(t), and P[q(t)] is the probability distribution of q at 
time t. The Bloch-Rcdficld equations, Eq.©, give the 
average time evolution of an ensemble of systems with 
the same initial condition, averaged over the details of 
the uncontrolled decoherence processes. Thus we find 
that (q(t))\ qj , is simply the solution of the Bloch-Redfield 
equations, given the initial value qj. More explicitly, we 
find the following expression: 
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{AN(t) - A [R(t) cos Of - I{t) sinOt]}, 

(11) 

where the initial condition qj is written in terms of the 
initial values Nq, i?o and Iq. The corresponding expres- 
sion for q(t) is 
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[AN SS - A (R ss cosM - I ss sin fit)] , (12) 



where N ss , R ss and I ss are the steady state limits for 
N(t), R(t) and I(t), respectively, obtained from Eq. ([7| 
at t — > oo. 

Thus we find that the correlator defined by Eq. (fTU)) 
depends on the phase of the driving field at both t\ and 
t 2 . In a qubit experiment, where t\ and t 2 are the initial- 
ization and measurement time, respectively, we assume 
that we do not have sufficient control over the phase of 
the driving field at initialization time and therefore aver- 
age over the initial phase of the driving field. The details 
of this procedure are given in Appendix [X] 

The details of the procedure used to calculate the two 
time correlation function, Eq. (fTDj) . by use of the Bloch- 
Rcdficld formalism is described in detail in Appendix IBl 
The full procedure includes a coordinate transform, and 
in the following we give an outline of the procedure. First 
we find the density matrix in the steady state that we 
might visualize as a point within the Bloch sphere. Next, 
we note that in an external field the steady state solution 
can in general lie anywhere in the Bloch sphere and not 
necessarily along the z-axis. Therefore, it is necessary 



to transform to the coordinate system where the steady 
state solution lies along the z-axis. The angles defining 
this transform are given by the steady state solution of 
Eqs.(6-8), and are illustrated in Fig. [TJ The angles and 
</> are defined by the relations: 
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The details of this transform and its application to the 
evaluation of the two time correlation function is de- 
scribed in Appendix [Bj Since the off-diagonal elements 
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FIG. 1: The coordinate transform used in order to diagonalize 
the density matrix. Here N, R and I are the parameters 
determining the density matrix in the energy eigenbasis of the 
TLS, while N' , R' and I' denote the same parameters in the 
rotated frame defined by Eqs. (|13[) . The frame is defined such 
that the off diagonal elements of the density matrix vanish 
in steady state, i.e., R' ss = /g S = 0. In external driving, the 
steady state values of the off-diagonal elements of the density 
matrix are in general non-zero. The transform is used in order 
to make use of the average procedure described in Appendix 
IBl In the absence of driving the two frames coincide. 

of the density matrix vanish in this choice of basis, we 
are allowed to use the states |— ') and |+') in the rotated 
basis as initial states, weighted by the mean population 
levels obtained from the density matrix in the steady 
state, that gives us P[q]. 

In the absence of external driving the density matrix 
in the equilibrium will always lie along the z-axis and the 
ensemble average of R and / vanishes. Thus we do not 
require the coordinate transform. In this particular case, 
after introducing r = t 2 — 1\, the explicit expression for 
the two-time correlation function given by Eq. (|10[) is 
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(14) 



Here An(r) denotes the 11 element of A(r) given by 
Eq. ((HJ). The dependence of N(t 2 ) on the initial values 
R(ti) and I(t\) vanish in the absence of external driving. 
We can therefore in this simple case write the propagator 
A(r) as a scalar function An(r). 

In the general case, when driving is included, we find 



G(h,t 2 ) <x f(T,A, 7 ,z)N^(l- N^ s ). 



(15) 
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Here /(r, A, 7, z) describes the dynamics of the density 
matrix, while N' ss is the population of the upper level in 
the rotated frame, illustrated in Fig.l. We have 

Ks^-Ks) = N sa (l-N sa )+g(A, 7> z) (N ss - 1/2) 2 (16) 

where g (A, 7, z) describes the details of steady state den- 
sity matrix (g(A, 7, z) = if R ss = I ss = 0). See Eq. (|B7j) 
and Eq. (|B8|) for details. 

From the correlation function, Eq. (|10|). we can com- 
pute the contribution of a single TLS to the noise spec- 
trum acting on the qubit. The spectrum is given by the 
expression: 



S(u) 



e^G{\T\,Q)dr, 



(17) 



where we took into account that the correlation func- 
tion G{t\,t2) is translation invariant after the averaging 
procedure described in Appendix [Bj We note that the ir- 
reversible Bloch-Rcdficld equations require the measure- 
ment time to succeed the preparation time, therefore we 
need the absolute value of |r| in the definition. From the 
full spectrum at arbitrary frequency, given by Eq. (|C1|) in 
Appendix [C] we obtain in the limit f2 > 7 the following 
expression for S(u>): 
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a 2 L(uj ± A') ± a 3 UJ±A ' L{u ± A') 



(18) 



For long measurement times this distribution becomes 
narrow, such that the pure dephasing rate of the qubit is 
given by = irS(0) when t — > 00. However, for re- 

alistic measurements with finite measurement times, the 
dephasing rate is determined by the noise spectrum in a 
finite domain of low frequencies centered at u> = 0. In 
the following, we will first derive expressions for the low- 
frequency contribution to S(uj) due to a single TLS in 
the absence of external driving (A —> 0) and in the case 
of strong driving (A » 7). These expressions will later 
be used to derive the low frequency noise spectrum for a 
specific distribution of driven TLSs. 

From Eq. (jl8p we find a crossover from a regime where 
the driving contributes as a weak perturbative effect to a 
regime strongly dependent on the driving, the crossover 
takes place around A w \z\, given that A ^> 7. We pro- 
ceed by deriving the limiting expressions in the resonant 
region \z\ <C A and in the off resonant region \z\ ^> A. 

Using the expression A oq = (e E l kT + l) and Eq. (fT8|) 
we get in the off resonant regime, \z\ A, 
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As follows from the above expression, the driving only 
weakly [oc (A/z) 2 <C 1] affects the fluctuators that are 
far from resonance. 

The corresponding leading contribution to the noise 
spectrum in the resonant regime \z\ < A, is 



where L(uj) = r )/{ r y 2 + w 2 ), 
oi 



z 2 cos 2 1 



zA sin 9 cos 9 cos (f>, 
d2 = A 2 cos 2 9 + zA sin 9 cos 9 sin (j>, 
03 = A A' sin 9 cos 9 cos <j>. 

In this limit, the spectrum only contains peaks at zero 
frequency uj = and at the Rabi frequency A' . 



B. Low frequency noise 

In the remaining part of this article we are interested in 
the noise at low frequency, and how it is changed by the 
external driving. The reason behind this focus is that 
the noise at low frequencies has been identified as the 
dominant source of pure dephasing in Josephson junc- 
tion qubits. It is known that for the Gaussian noise and 
diagonal qubit fluctuator coupling, the off-diagonal ele- 
ments of the qubit density matrix relaxes (in the case of 
the free induction decay) at a rate (see, e.g., Ref. [l2t) 
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From the full spectrum Eq. (|C1|) . together with Eq. (jT6|) 
for the population of the density matrix, we can identify 
two main effects of external driving on the TLS noise 
spectra. The first effect is the altered equilibrium pop- 
ulation of the density matrix due to driving. We can, 
by use of Eq. ©, express the occupation of the density 
matrix in steady state by 



N, 
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1 + (A 2 / 7 )L(2 
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as previously found in Ref. [Tj. 

External driving results in saturation of the steady 
state density matrix when A sa 7 for \z\ <C 7 and when 
A w \z\ for \z\ ^> 7. This saturation contributes to 
increased fluctuation rate between the upper and lower 
level of the TLS, see Fig. [51 and therefore this effect con- 
tributes to a higher intensity of the noise at all frequencies 
(that was not found in Ref. [HI). This is especially true 
for TLSs where the energy splitting is large compared 
to temperature, meaning that the noise is very weak in 
thermal equilibrium since the system spends almost all its 
time in the ground state (see Fig. [3]). There is, however, 
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FIG. 2: Noise spectrum induced by a single TLS for A — 
(solid line), 5 x 10 6 Hz (dashed line) and 5 x 10 7 Hz (dash- 
dotted line) . We see that with this choise of parameters, driv- 
ing reduces the noise at zero frequency, but enhances the noise 
at higher frequencies. The peak at the renormalized Rabi fre- 
quency, A' , is the most pronounced. Peaks at higher frequen- 
cies are suppressed as long as A <C f2. The parameters used in 
the figure are T = 0.2 K, S7 = E/h = 10 10 Hz, 7 = 5xl0 6 s _1 , 
70 = 10 7 s" 1 , 
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FIG. 3: Noise spectrum induced by a single TLS, with the 
same parameters as used in Fig[2] but at lower temperature 
T = 0.05 K. At this temperature, the noise is weak in the 
absence of the driving (solid line) since the fluctuator is frozen 
in its ground state, N cq ~ 0.01. Here the main effect of 
driving with strenght A — 5 x 10 6 Hz (dashed) is to increase 
the probabiliy of excitation leading to increased noise at all 
frequencies. When the strength of the driving is A — 5 x 10 7 
Hz (dashed/dots) the noise spectrum is shifted sufficiently 
away from u> — towards cj = A' such that the effect of 
increased fluctuations is offset by the shift in the spectrum. 
Thus the noise at u = is reduced. 



another very pronounced effect, not caught by the model 
of Ref. IH. The full noise spectrum of the driven TLS 
Eq. (jCip is composed of several peaks. At low frequencies 
the most pronounced are the one centered around uj = 
and the two peaks centered around the renormalized Rabi 
frequency, uj = ±A', see Fig. [5J From the last term of 
Eq. (fT5)) we see that the intensity around the uj rs peak 
is reduced when the Rabi frequency A becomes compa- 
rable in magnitude to the deviation from resonance z. 
Thus we have a shift in the intensity from low frequen- 
cies to frequencies around the Rabi frequency. This shift 
might be beneficial in reducing the low frequency noise 



responsible for pure dephasing of the qubit In a sim- 
plified picture the total noise at a given frequency as a 
response to external driving can therefore be regarded as 
a result of two competing mechanisms, the increased fluc- 
tuations due to increased population in the upper level of 
the TLS, and the shift in the spectrum from low to high 
frequencies. At sufficiently strong driving this shift leads 
to a reduction in the low-frequency (w < 7) contribution 
to the noise spectrum from a TLS with E < kT as 



S(u < 7) 



(res) 
A>7 



S(u < 7)0 



57 2 
2A* 



(23) 



IV. ENSEMBLE OF TLSs 

In this section we will analyze the noise from an en- 
semble of the fluctuators studied in the preceding sec- 
tion. Our purpose is to roughly estimate the effect 
of external driving on the noise spectrum for a realis- 
tic Josephson-junction qubit experiment. We will more 
specifically assume the following distribution of the TLS 
parameters. First, we note again that our calculations 
are based on the assumption that the dynamics of the 
fluctuators are independent of the state of the qubit. Fur- 
thermore, we assume in the following that the strength 
of the fluctuator-qubit coupling v is uncorrelated with 
the relaxation rate 7 and the energy E. Assuming that 
Ao is an exponential function of an almost uniformly dis- 
tributed parameter, such as tunnel barrier height ; 12 i 26 i 27 



the distribution of the TLS parameters becomes 
P(A,A ) = Ptls/A 



(24) 



where Ptls is proportional to the density of states per 
unit energy and volume. This distribution is already 
widely used in models of decoherence in qubits, where 
it is known to give origin to the oc 1/cj dependence of the 
noise spectrum at low frequencies. 

Since in the following it is more convenient to work 
with the relaxation rate 7 and the unperturbed fluctua- 
tor energy E, we recast Eq. (|2"4"]l by use of the relationship 
7 = 7o(P)(A /P) 2 , where jo(E) is the maximum relax- 
ation rate for a fluctuator of energy E, obtaining 



Ptls 



7\/! - 7/7o(P) 



(25) 



for 7 <E [7min, 7o(P)]- The distribution has to be cut 
at the relaxation rate 7 m i n of the slowest fluctuator. 
However, we find that the noise spectra at frequencies 
uj 3> 7min; and therefore measurements carried out with 
the measurement time r = 1/ui <C l/7min, are not sensi- 
tive to the cutoff. The maximal relaxation rate, 70, is a 
power-law function of the energy E. Since in the follow- 
ing we restrict ourselves to order-of-magnitude estimates 
we will replace 70 (E) by a constant rate, 70 ~ 70 (&T). 
Our calculations (see Appendix [TJ]) show that using this 



7 



assumption the noise spectrum depends only weakly on 
To- 

Before we proceed to evaluation of ensemble integrals 
it is convenient to introduce a new variable 



= E ac p/ VTo- 



(26) 



This variable is independent of 7 and can be treated like 
a constant when integrating over distributions of TLSs. 

If we assume that the qubit-fluctuator coupling v is 
uncorrelated with 7 and E, and by using that the single 
fluctuator spectrum is cx v 2 , we can express the ensemble 
averaged noise spectrum by 



70 



S a (u) = (V 2 ) 



S a (uj,E, 7 )P(E, 7 )d 1 dE. (27) 



7„ 



Here we have introduced the notation (v 2 ) = 

Umax 

J v 2 P(v)dv, where P(v) is the distribution of the 
qubit-fluctuator coupling v. By use of the given distribu- 
Y 




FIG. 4: The full TLS parameter domain. The TLSs inside 
the parabola (fi — E/K) 2 < A 2 , group II, are resonant with 
the applied field. For driving frequencies Q > kT/h the major 
contribution to noise origin from low-i? fluctuators belonging 
to group I outside the resonant sector. The contribution to 
the noise due to these fluctuators are not changed significantly 
by driving at frequencies much higher than their energy split- 
tings E. 



tion of TLSs, Eq. (|25j) . and the expression for the noise 
in zero driving, Eq. (|20[) , we can evaluate the ensemble 
averaged noise in the case of no driving (a = 0). The 
detailed calculation is given in Appendix [DJ Using the 
relationship A 2 /E 2 = 1 — 7/70 between the relaxation 
rates and the fluctuator potential parameters^ we find 
that the averaged spectral density is given by 



So(w) 



1 u \ 7 min < u < 70, , , 
v )kTP TLS i _ -1 , . <r- ( 28 ) 



Tn 



w i$ T 



We conclude that without driving we obtain noise cx 1/uj 
for the interval 7 m i n < uj < 70 that turns over to a con- 
stant value for w < 7 m ; n . 



Next we proceed to strong driving, which we have de- 
fined by a 2 > 7o- In order to evaluate the noise in this 
regime we split the domain of integration in two parts, see 
Fig. 2] The resonant domain a 2 j > z 2 (group II), where 
the fluctuators are strongly affected by the external field 
E ac , and the offresonant domain a 2 j < z 2 (group /), 
where the fluctuators does only weakly respond to the 
driving field. We approximate the full integral by us- 
ing the asymptotic limits given by the undriven noise 
spectra, Eq. ([2U]). in the offresonant domain, while the 
strong driving limit is given by Eq. (|2ip . The total en- 
semble averaged noise spectra for frequencies in the in- 
terval 7 m i n < uj < 70 is (see Appendix [D] for details of 
derivation) 



-{v z )Ptls 

7T 



UJ 



Aha 



1 



cosh 



2 m_ 

2kT 



(29) 



for hciy/uj <C kT and y/jo < a < fl/y/jo. Again, at 
uj < 7min the frequency uj in this expression should be 
replaced by 7 min . 

From this result, we find a correction to the 1/uj noise. 
This correction, being oc uj^ 1 / 2 , is due to suppression 
of the noise from TLSs that are close to resonance with 
the driving field f2. The driving therefore leads to a re- 
duction in the noise at intermediate and low frequen- 
cies. The derivation of the expressions Eqs. (|2"9"]l and a 
similar for uj < 7min is a good approximation only as 
long as a <C ksT/hy/jmin. As long as this criterium 
holds, we find that the correction term is small relative 
to the first term, even if the ratio KVt/kT is small and 
the term cosh 2 (ft£l/2fcT) approaches unity. If, however, 
we increase a beyond this inequality, we expect that the 
correction term will increase until it approaches the first 
term in magnitude. The physics is as follows. When we 
increase a we increase the number of fluctuators belong- 
ing to group iT, at the cost of group / by increasing the 
width of the parabola z 2 — A 2 in Fig. @] until all fluc- 
tuators are resonant with the field. Thus the number 
of fluctuators responding to the driving field is increased 
and since each fluctuator within the resonant sector will 
have its noise spectrum shifted toward higher frequen- 
cies, the noise at low frequencies will be reduced until it 
approaches zero for very large fields. 



V. DECOHERENCE OF THE QUBIT 

In this section we will describe the decoherence of the 
qubit due to the driven TLSs in its environment. We will 
illustrate the effect from driven TLSs by use of an exam- 
ple with a specific, but motivated distribution of TLS 
parameters. We will still assume that the qubit couples 
diagonally to the TLSs. The coupling was previously 
specified to be Hp_q = VifX z eg) t z . Given this coupling 
the TLSs does only have a pure dephasing effect due to 
renormalization of the qubit level splitting, and direct 
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transitions between the levels of the qubit (Ti processes) 
cannot be induced by our two level fluctuators. 

Above we have assumed that the TLSs-qubit coupling 
Vi is not correlated with 7 and E. However one have 
to keep in mind that different distributions of Vi might 
have significant impact upon the dephasing of the qubit. 
For in-depth treatment of different ensembles of fluctu- 
ators, as well as non- Gaussian noise statistics, we refer 
to Refs. [lO. 12. 2 8. 2 ij. Since we are in this article primar- 
ily interested in the effect of external driving, we assume 
that all fluctuators couple to the qubit with the same 
strength v. 

In the standard Gaussian approximation, the pure de- 
phasing time T2 for long times t is approximately given 
by T 2 T 1 = nS(Q), where S(0) is the noise spectrum at zero 
frequency^ Using this formula, together with Eq. (|D10|) 
for the noise at frequencies u < 7min, we obtain the fol- 
lowing expression for the dephasing time of the qubit 

1 , — , „ / kpiT Aha \ 

T 2 \^7mm ^7~cosh J 

valid for fi.a7 m i n <C kT and y^To < a < Q/ y/jo. 

From this expression we see that if Ml > kT, then 
the function cosh 2 (fifi/2jfcT) sa fW/KT ^ Wc thcn find 
that the relative reduction in the noise spectrum due to 
driving is oc ha^/j m i n / kTe hn / kT . Thus the correction is 
exponentially suppressed at low temperatures. However 
we notice that if the driving frequency is reduced, i.e., if 
one introduce a driving field at a frequency Ml < ksT 
much lower than the qubit frequency, the correction due 
to driving will become significant and the noise at zero 
frequency will be reduced. 

VI. DISCUSSION 

Our main result in this article, is that for Joseph- 
son junction qubits where the dominant noise source is 
TLSs interacting with the qubit, external driving have 
two main effects. The first is increased fluctuations of 
the TLS, contributing to increased noise at all frequen- 
cies. This effect is significant if the energy splitting of the 
fluctuator is small compared to temperature. The second 
effect is a reduction in the noise spectrum at low frequen- 
cies and at the same time increased noise at high frequen- 
cies, especially at the renormalized Rabi frequency A' . 

For a typical ensemble of fluctuators characterized by 
the distribution P(A, Ao) ex 1/Ao we find that external 
driving at high frequencies (e.g., the qubit frequency) 
leads to reduced noise spectra at low to intermediate fre- 
quencies, which again result in an enhanced qubit de- 
phasing time 2~2 • For a typical distribution of TLSs, 
P(E,j) oc I/7, the effect is weak since the driving only 
reduces the noise from TLSs with energy splittings close 
to resonance with the driving field (group II). For driv- 
ing fields Ml > kT, the resonant fluctuators (group II) 
contribute only weakly to the noise also in the absence of 



driving. In this case, both the driven and the undrivcn 
noise spectra at low frequencies is strongly dominated 
by TLSs (group /) with low energy splittings E < kT 
and long relaxation times T\ = 1/7 < uj^ 1 . However, 
by driving at a lower frequency Ml < kT, the resonant 
group (group II) is shifted from fluctuators close to the 
qubit frequency to fluctuators at lower frequencies that 
contribute stronger to the dephasing noise on the qubit. 
In this case we expect a strong reduction of the low fre- 
quency noise. 

An important side effect of external driving of the TLSs 
that we can identify in Eq. (|C1[) is increased high fre- 
quency noise around the renormalized Rabi frequency A' , 
the driving frequency fl and at beatings between these 
two frequencies. In experiments where the relaxation 
time is an important limiting factor in preserving 
the coherence of the qubit we expect that the increased 
high-frequency noise due to external driving will be coun- 
terproductive and one should take measures in order to 
make sure that the high frequency peaks does not overlap 
with the qubit frequency. 

Appendix A: Averaging over the phase of the 
driving field 

The expression for the correlation function in the time 
domain, Eq. (|10[) , can be written out explicitly by use of 
Eqs. (JTTJ and flU}. We obtain: 

G(h,t 2 ) = (F(t 2 )F(h)) , (Al) 
F(t) = ASN(t)-A [SR(t) cos fit - SI(t) sin fit] . 

Here we have introduced the notations SN(t) = N(t) — 
N ss , SR(t) = R(t) - R ss and SI(t) = I(t) - I ss . This 
expression depends explicitly on the exact phase of the 
driving field, E ac , at times t 2 and t\. In the following wc 
assume that the phase of the field is random at the start 
of the pulse. This is typically the case if the rise time 
of the signal is long compared to the oscillation period 
27r/0 of the signal. In order to average over repeated ex- 
periments with a random distribution of the phase of the 
field, we replace flt\ by flt\ + ft, where we assume that ft 
is uniformly distributed on the interval ft G [0, 27r]. With 
this assumption only products of sines and cosines con- 
tribute, while single terms vanish after averaging over the 
phase ft. After this averaging procedure we find that the 
expression for the two time correlation function Eq. (|A2[) 
is reduced to 

G(t, f + r) = (4v 2 /E 2 ) (A 2 6N(t + r)SN(t) 

+ (A 2 /2) {SR(t + t) [5R(t) cos Sir + SI(t) sin Sir] 
+81 (t + t) [SR(t) sin Qt - SI (t) cos Sir}}) (A2) 

where the cross terms of Eq. (|A2[) proportional to AAo 
have canceled due to averaging over repeated experi- 
ments with random distribution of the initial phase of 
the driving field. 
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Appendix B: Coordinate transformation and 
averaging over initial conditions 

The Bloch-Redfield equations are equations of motion 
for the average of an ensemble of TLSs, where the indi- 
vidual members of the ensemble differ by details of the 
environment. By use of the Bloch-Redfield equations we 
avoid dealing with these details that we do not have con- 
trol over. In place we get an equation of motion for the 
mean density matrix. While this method greatly sim- 
plifies the dynamics, since we are no longer required to 
keep track of fine details of the environment, the cost is 
loss of information about the time evolution of individual 
systems. Therefore the Bloch-Redfield equations cannot 
be used to calculate two time correlation functions in a 
straightforward way 30 ' 31 The procedure we use in order 
to evaluate the two-time correlation function is the fol- 
lowing. In general the two-time correlation function can 
be expressed as 

(i4(ti)A(t 2 ))=53o*(t a )C(o*(*2)|a J (t 1 ))oj(*i)C(ai(t l )). 
j,k 

Here A(t) and a(t) is an observable and a particular re- 
alization of this observable, respectively, C(ak(t2)\ aj (t 1 )) 



is the conditional probability distribution for observ- 
ing the value afc(£ 2 ) at time i 2 conditioned upon that 
the value Oj(ti) was observed at t\. And the £(aj(ii)) 
is simply the probability distribution for observing the 
value aj(ti) at time t\. We might then realize that 
E a fc( t 2)C( a fc(i2)U 3 (t 1 )) = (^(*2))L(ti) is simply the so- 

k 

lution of the Bloch-Redfield equation given the initial 
value aj{t\). The two-time correlation function thus re- 
duces to 



{A(ti)A(t 2 )) = ^(A(t 2 ))|^ (tl) a J -(t 1 )^(a i (t 1 )). (Bl) 
j 



If we now move to our specific problem of a driven TLS 
in a dissipativc environment, the two time correlation 
function we require, Eq. (jlOp . does only contain terms 
with the deviation of the observable quantity from its 
steady state value. From the explicit solutions of the 
Bloch-Redfield equations, Eq. (|8]) and Eq. (|9]), we see 
that the time evolution of the deviation from the steady 
state is translation invariant, linear and homogeneous. 
Therefore, we can write 



(9(*a) - 9(*2)>| w (t0-5(tx) - (2d/E){[AX n (r) - A (A 21 (r) cos fir - A 3 i(r) sin fir)] (A, (t^ - N ss ) + [AA 12 (r) (B2) 
-A (A 22 (r) cos fir - A 32 (r) sinfir)](i^(ti) - R ss ) + [AAi 3 (r) - A (A 23 (r) cos fir - A 33 (r) sin fir)] (Ij (t x ) - 7 S8 )}. 

I 



Here r = i 2 — 1\ and A Q( g(r) arc the elements of A(r) given 
by Eq ([5]). This expression can be inserted directly into 
our general formula given by Eq. (|B1|) . Unfortunately 
we still require the distribution function £(Nj(ti)), and 
similarly, the distribution of / and R. The distribution 
functions cannot be extracted from the Bloch-Redfield 
equations^. To make up for our lack of information, 
we make the following approximation. Assuming that a 
measurement of the TLS in the eigenbasis of i?o, will give 
either the outcome Nj—o = or Nj—i = 1 with the mean 
value A ss , we find that 

Njt{Nj) = n q z(No) + mam) = A ss , 

3 

]T Nft(N,) = Nfe(No) + N^im) = N ss . (B3) 
j 

In addition, with this choice of initial values, the initial 
values of the off-diagonal density matrix elements is al- 
ways zero. We note that this choice of initial values only 
make sense if the steady state density matrix lie on the 
axis between the points Nj—o — or Nj—i = 1. In order 
to apply the method, we are therefore required to trans- 
form to the coordinate system where the steady state 



value of p lie on the z-axis of the Bloch sphere. 



1. Coordinate transformation 

The summation procedure derived above, Eq. (|B3I) . 
works nicely as long as the assumptions behind the 
derivation of the Bloch-Redfield equations (the Born- 
Markov approximations) are fulfille d 20 i 21 , as well as the 
time translation invariance. However, the summation 
procedure can only be applied if the steady state of the 
density matrix lies along the z-axis of the relevant mea- 
surement operator. With driving, we see that the density 
matrix will in general be driven away from the z-axis, 
such that a summation over the eigenstates |+) and |— ) 
of the Ho operator given by Eq. ([2]) , can not possibly give 
the true average density matrix in steady state. However, 
we can do proper averaging by transforming to a new co- 
ordinate system where N' ++ denotes the occupation along 
the z'-axis in this new choice of coordinates. After this 
transformation the Bloch vector of the steady state den- 
sity matrix N ss is a point on this axis. By this choisc 
of axes the density matrix is diagonal. The coordinate 



10 



transformation is given as 

(N - 1/2) = (N' - 1/2) cos6, 

R = N' sin 9 sin <p, I = N' sin 9 cos <j>. (B4) 

We note that R' sa = and I s ' s = since we have defined 
Ng S to lie on the z'-axis in the new coordinate system. 
When we insert the explicit steady state expressions into 



Eq. (JT3J), the dependence on the equilibrium value N eq 
vanishes, and the expressions reduce to 



tan# = 



q 2 7 (z 2 + 7 2) 
4( 7 2 + z 2 + 2a 2 7 ) 2 ' 



tan </> = — 

7 



(B5) 



Inserted into our expression for the correlation function 
Eq. (|A2[) we obtain the following formula: 



G{hM) = (2v 2 /E 2 )(2A 2 6N(t 1 )[X 11 (T)5N(t 1 ) + Ai 2 (r)Mi(fi) + Ai 3 (r)*J(ti)] + Ag{[A 2 i(T)<H\r(ti) + A 2a (T)MJ(ti) 
+A 23 (T)<fJ(fi)][Mi(ii)cosfiT + <5I(ti) sin Or] + [A 3 i(T)(JiV(ti) + X 32 {T)5R(ti) + A 33 (r)W(ti)] 

x cos Or + sin Ht] }) . (B6) 



We can now by use of Eq. (|B4[) move to the frame where the density matrix lie along the z-axis. By use of the 
summation formulas given by Eq. (|B1[) and Eq. (|B3[) . and inserting the explicit expressions for the elements of A(t), 
given by Eq. ([8]) we obtain the following expression for the correlation function: 

G(t, 0) = (2v/EA'f iV4(l - 7V S ' S ) { A 2 [(A 2 cos A't + z 2 ) cos 2 9 - zA{l - cos A't) cos 6 sin 6> sin 

A 2 

+AA' sin A't cos 6* sin cos 0] H — — [A' 2 sin 2 9 cos cos A't cos Qt (cos — sin cf>) + AA' cos 9 sin 6> sin A't 

z A 2 
x (sin sin fir — cos cos fir) + — sin 2 9 sin A't sin Ht + — ^ sin 2 sin 0(cos sin J7t — sin </> cos Ot) 

+z sin 2 6* sin <p sin A't(cos </> sin fir — sin </> cos Qt) — zA cos 9 sin 6*(cos </> sin fir + sin cos Ot) 

+zA cos 6* sin cos A't(cos </> sin fir — sin cos Ot)] } . (B7) 

Here 

<(1 - O = ^ss(l - N ss ) + (n ss - sin 2 - ^ sin 2 9 sin 2 - ^ sin 2 9 cos 2 

1\~ , ( , - 1\ T , „ , . „ , ^ss *^ss . 2 



iV ss — — J R ss cos 6* sin sin (f> — yN ss — — J I ss cos(#) sin 9 cos <f> — sin 9 cos <fi sin <f>. (B8) 

Eq. (|B7[) is our final expression for the two time correlation function for a single nuctuator. We note that the 
correlation function has become fully translation invariant after the averaging procedure. The Fourier spectrum is 
computed in Appendix [Cl 



Appendix C: Spectral density S(cj) at arbitrary frequency 

In this appendix we give the full spectral density from a single TLS. The full spectrum is given by the fouricr 

_ 00 

transform with respect to t = i 2 — t\ of Eq. (|B7[) . Thus we calculate S(lo) = J e lUT G(\r\, 0)dr. Carrying out 

— oo 
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the transform wc obtain the result 



S(oj) =8^ (^7) 2 K a (l - iV s ' s )(A 2 {A 2 cos 2 9 + z A sin9 cos 6 sin 0[L(w + A') + L(oj - A')] 



~AA' sin 9 cos 9 cos < 



A' + ui 



A' 



L(u) + A') + — -L(u -A') + (z 2 cos 2 9 - zA sin 9 cos 9 sin 

7 7 J 



A 2 r 

+ ^ sin 0|&j_ [L(w + A' + fi) + L(u - A' - fi)] + b 1+ [L(lu + A' - fi) + L(w - A' + 0)] 



-&2+ 
-b 2 - 



7 7 

A'-n + uj. n . -A' + Q, + lj t . ., 
L(w + A — fi) H L(u) -A' + Sl) 



+2b 3 j sin^[L(w + fi) + L(w - fi)] + cos0 



fi 



±^L(o; + fi) + ^L( W -fi)]}}), 



(CI) 



A' 2 sin 9 cos 0(cos — sin (j)) ± AA' cos 9(1 + cos<p) + z sin 9 + z 2 sin 9 sin + Az cos 0) (cos + sin < 
—AA' cos 9 cos 4>± z 2 sin 9 sin </> ± Az cos 0) (cos <fi + sin <^>) , 
A 2 sin sin </> — Az cos 9. 



(C2) 



From Eq. (|C1[) . we see that without driving (and the 
same for off- resonant driving), we have a single peak 



S(ui) oc 



1 A sin 6* cos 9 cos <j> 7 



A 2 + z 2 



7 2 + w 2 



that is reduced as the driving increases (i.e., when the 
Rabi frequency A approaches z). When the driving is 
strong, the intensity is shifted from the single peak at 
lu = 0, to a large number of peaks. Most prominent are 
the peaks at the renormalized Rabi frequency, uj = A', 
but there are also peaks at the driving frequency, and at 
sums and differences between A' and fi. The noise at 
bj « fi has not been discussed in this article, but might 
be important for qubit operation if the driving field is 
the manipulating pulses used to control the qubit. This 
noise is then close to resonance with the qubit. 



spectrum in resonant strong driving given by Eq. (|2 1 [) . 
while in the off resonant sector we use the expression for 
fluctuators out of resonance given by (|2"0|) . Our strat- 



Yo 



z 2 =A 2 




FIG. 5: The domains of integration 



Appendix D: Derivation of the ensemble averaged 
noise for strong driving 

In this section we will derive the ensemble averaged 
spectrum of the low frequency noise induced by strongly 
driven TLSs (i.e., a 2 > 70) for the distribution of TLS 
parameters given by Eq. (|25[) . In order to evaluate the 
integral over the TLS parameters we make the follow- 
ing approximation. The full parameter domain 7 ® E £ 
[7min, 7o] ® [0, 00] is split in two sectors. In the first sec- 
tor the TLSs are in resonance with the driving field de- 
fined by the criteria a 2 7 > z 2 (group II). The second 
one contains the TLSs that are out of resonance (group 
J, defined by a 2 7 < z 2 ). The domains of integration 
is specified in Fig. [5J For the fluctuators belonging to 
the resonant sector we use the expression for the noise 



egy to evaluate the noise spectra is to first compute the 
ensemble averaged noise spectrum in the absence of driv- 
ing over the full domain, then subtract the contribution 
to the undriven spectrum from fluctuators that lie in the 
resonant sector, and finally we add the contribution from 
the resonant fluctuators in group II. Therefore, the noise 
spectrum can be represented as 



S, 



j(res) 



(Dl) 



Here the superscript denote the resonant sector. 

We start by evaluating the spectrum in the absence of 
driving, by use of Eqs. (j20|) and (f25|) we get 



5oH = A /-V) 



00 70 



1-^- 

70 



7 



TLS V 

, 2 - 2 , r 2 e dEd ^ (° 2 ) 

(7 2 +cj 2 )cosh 
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For 7min *C oj *C 70 the integral over 7 can be estimated 
as 



1 , 7o 
■In — , 



T 



uj ui 2 70 270 oj 



where we have used the expression for the resonant spec- 
trum in strong driving given by Eq. (|2ip . This integral 
is evaluated similarly to the corresponding integral for 
the same region in the absence of driving. After using 
asymptotic expressions in the different regions of Fig. [5j 
we find 



while the energy integral is J dE cosh 2 (<d|f) = 2kT. 



Therefore 

S (oj) 

where A = \J 
w < 7min we get 



M 2 - 

■.UJ 



1 

7o 



1 

2^ 



In 2° 



1 'n '(v 2 ) Ptls^T ■ For very low frequencies 



S a=0 {uj) ~ A 



hi 



1 

To 



-!-ln — 

270 7min 



(D3) 



Next we proceed to calculate the contribution from the 
resonant sector in the absence of driving. The integral 
we need to evaluate is similar to that of Eq. (|D2[) . but 
the integration should be performed over the TLSs from 
group II. The integration domain is restricted by the 
parabola z 2 < a 2 j. 

Assuming first that 7 m j n < oj < 70, we evaluate the 
integral by use of asymptotic expressions for the inte- 
grand in the different regions of Fig. [5] More precisely 



we use 



-j2 in sector A of FigJSJ and 



1 

7 2 



in sector B. Furthermore we make the approximation 



cosh 



2 e 



cosh 



2 HQ 



- const for all fluctuators in- 
side the resonant sector. If we write E = h(il — z), we 
find that since z 2 < a 2 70 inside the resonant sector, the 



approximation is good as long as -^j, < jj T <C 1. 
However, the major contribution to the integral comes 
from 7 < oj, such that we can narrow our inequality to 



-np _^ kT -C 1. By use of the approximations described 
above, we find that the total contribution from the res- 
onant sector in the absence of driving, in the frequency 
interval 7 m j n < w < 70, is 



S' 1 



(res) 



A 



2ha 



kT cosh 2 ^\ 




For uj < 7 m in, the calculations are similar, but somewhat 
simpler. The result is 



S, 



(res) 



A- 



2ha 



V7n 



kT cosh 2 V Vluln' y/jO 
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(D5) 

Finally we proceed to the resonant sector in strong 
external driving, a 2 > 7: 



Ares) 
V> 7 



5.4 

2a 2 kT 



1 



7o 



djdE, 



(D6) 



group II 



Ares) 
'a 2 > 7 



25Ah 



WakT 



(D7) 



for the total contribution from the resonant sector in ex- 
ternal field. This expression is valid for oj < 70- 

We have now computed the three contributions to 
the total ensemble averaged noise spectrum given by 
Eq. (|Dip . For 7 m ; n < oj < 70 the full spectrum is given 
by 



J a*>-y 



A 



/_2 _ 7mm _ _1 Lin^^l 

\u> oj 2 70 270 uj ) 
2ha 



CO OJ 

25V70 
4afcT 

3 



kT cosh 2 Jjg, \3y/UJ 



2 7 3/2 

[ mm 

3w 2 



7o 



(D8) 



Within the limits of the inequality used for the evaluation 
of the integrals, hay/oo <C ksT, and by using that a > 
y/ui, we see that the second term origining from the driven 
resonant sector is negligible compared to the undriven 
1 /oj term. Thus we are left with 



2ha 


8 


kTcosh 2 ^ 


\3V^ 




For frequencies uj < 7 m in the contributions to the en- 
semble averaged spectrum is given by Eqs. (lD3p , (|D5|) and 
(|D7]) . giving: 



5, 



2ha 



A 



1 



kT cosh 



2 m 

2kT 



7m in 



70 27o 7min^ 



70 



(D10) 
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